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^ We explore the dark matter interpretation of the anomahes claimed by the DAMA 

,__! and CoGeNT experiments, in conjunction with the various null direct-detection ex- 

P^ periments. An independent analysis of the CoGeNT data is employed and several 

1^ experimental and astrophysical uncertainties are considered. Various phenomeno- 

^ logical models are studied, including isospin violating interactions, momentum- 

^— I dependent form factors, velocity-dependent form factors, inelastic scatterings (en- 

_ r*" dothermic and exothermic) and channeling. We find that the severe tension between 

/\ the anomalies and the null results can be ameliorated but not eliminated, unless 
extreme assumptions are made. 
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1 Introduction 

Significant on-going efforts are being made to directly search for Weakly Interacting Massive 
Particles (WIMPs). An interesting and widely accepted signature of WlMPs is the annual 
modulation of their interaction rate, arising from the relative motion of the Earth around the 
Sun. Strikingly, both the DAMA [1] and CoGeNT [2, 3] collaborations observe anomalous 
modulating events that may be interpreted as arising from interactions of spin-independent 
Dark Matter (DM). As we discuss below, the two measurements are, to some extent, consistent 
with each other, and point to a surprisingly low DM mass, of order a few GeV. 

In contrast to the positive signals of DAMA and CoGeNT, several other experiments 
find no evidence for DM. Most notably, the CDMS [4], XenonIO [5] and the XenonIOO [6] 



collaborations seem to disfavor the parameter space indicated by DAMA and CoGeNT. It is 
natural to ask, therefore, what possible systematic effects and/or DM properties can resolve 
the tension between the positive and null results? 

One noticeable difference between DAMA/CoGeNT and the null experiments, is that the 
latter veto electronic interactions while attempting to collect only nuclear recoil events. It 
is conceivable that the anomalous signals arise from such electronic recoils, a possibility that 
would explain away the existing tension. A model of this type was considered in [7, 8] prior 
to the recent CoGeNT measurement [3] and it remains to be seen whether this possibility 
is theoretically feasible. In this paper we pursue a different direction and study the viability 
of nuclear recoils of spin independent DM, as an explanation to the positive signals. Several 
handles can, in principle, ameliorate the tension with the null results. From the DM perspective 
one may consider, 

• Inelastic scattering (endothermic or exothermic) [9, 10, 11]. 

• Isospin-violating couplings [12, 13, 14, 15, 16, 17, 18]. 

• Velocity suppressed interactions. 

• Momentum dependent scattering [19, 20]. 

• Resonant scattering [21]. 

In addition there are uncertainties that may significantly change the expected scattering rates 
in various experiments: 

• Astrophysical uncertainties [22, 23, 24, 25, 26, 27, 28, 29, 30]. 

• Possible channeling effects [31, 32] (see however [33]). 

• Uncertainties in the Xenon scintillation function at low recoil energies or in the statistical 
treatment of the background events (for a discussion, see e.g. [34]). 

A fully systematic analysis that takes into account all of the theoretical and experimental 
uncertainties above is hard to attain and we do not attempt here. Instead, we separately 
study the effects of most of the above possibilities on the DM interpretation of the DAMA and 
CoGeNT results, as well as the null experiments. 

Our goal is two-fold. Primarily, we aim at understanding to what extent the DAMA and 
CoGeNT results are consistent with each other, and with the other null experiments (for 
related works, see [35]). In addition, we study what is required from the theoretical point of 
view and what needs to be assumed on the experimental side, in order to ameliorate the tension. 
Along the way, we reanalyze the CoGeNT modulation data. 

The paper is organized as follows. In Section 2 we summarize the standard computation 
of the expected signals and shortly discuss the DM velocity distributions considered here. In 
Section 3 we describe the experimental data we use, dwelling on the various uncertainties and 



discussing their influence on the fits. Our results are presented in Section 4. Here we follow 
the analyses endorsed by the various experimental collaborations and explore whether velocity 
distributions or any of the DM scenarios mentioned above can explain the signals compatibly 
with the bounds. We conclude in section 5. 

2 Formalism 

The direct detection rate for DM-nucleus scattering at a given experiment is given by 

dR ^r Po) f ,-i n , ^ da ,„ ^, 

-Nt^ dhvUv,t)-—, (2.1) 



dE^, Mdm J\v\>v^i,^ dE-R 

where Nt is the number of target nuclei per unit mass of the detector, Mdm is the DM mass, 
Pq is the local DM density (that we assume to be equal to 0.3 GeV/cm^) and fi^{v,t) is the 
DM velocity distribution in the Earth frame, to be discussed below. We denote a to be the 
DM-nucleus scattering cross section, which we take to be of the form, 

da _ rUNan [fpZ + fnjA - Z)f ^ . x ^2 . x .9 9^ 

Here m^ is the nucleus mass, /i„ is the DM-nucleon reduced mass, fp (fn) is the coupling 
strength to the proton (neutron) and F]sf{q) is the nucleus form factor. Throughout this work 
we use the Helm form factor given in [36]. The DM form factor, -Fdm(q'; ^); is a velocity and/or 
momentum dependent contribution to the cross section which may exist, depending on the 
DM coupling. Below we consider several possibilities for its form. Finally, fmin appearing in 
Eq. (2.1) is the minimal DM velocity needed for a scattering with a recoil energy E-^ to occur. 
In the elastic scattering case, it is given by 



'^. (2,3) 

In Section 4.8 we consider inelastic scattering, in which case the required minimal velocity 
depends further on the mass sphtting, 6 = M^j^ — Mdm, of the recoiling particles: 
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The scattering rate is expected to be time-dependent, exhibiting a maximum at the begin- 
ning of June (when the Earth moves against the DM wind) and a minimum at the beginning of 
December (when the Earth moves along the DM wind). It is thus useful to define a modulated 

rate as 

WR , 1 r WR WR 1 

(2.5) 



dET, 2 



dR , -, ^ dR / _^ , ^ 

(2 June) — -77:;- (2 December) 
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Both the modulated and unmodulated rates can be confronted with the results of an exper- 
iment only after accounting for all expected backgrounds and experimental efficiencies. The 
modulated rate is, in particular, independent of any time-independent background which can 
contaminate the DM signal. 
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Figure 1: Left: Dark Matter velocity distributions: the Maxwell-Boltzmann with sharp cutoff 
at V = fesc = 500km/s (thick black curve), the distribution, Eq. (2.8), with a smooth cutoff 
computed for k = 1 (red) and k = 3 (blue). Two different values of vq are shown: 220 km/s 
(solid) and 270 km/s (dotted). Right: The function rjn defined in Eq. (2.9) as a function of 
the minimal DM velocity, fmin; normalized to the mean DM velocity, vq. As discussed in the 
text, this function captures the effect of velocity- dependent form factor on the total scattering 
rate. We plot F^^{v) = 1,v'^/vq,v^/vq (magenta, blue and green respectively) for the Maxwell- 
Boltzmann distribution. For completeness we also show the corresponding function for the fk 
velocity distribution defined in Eq. (2.8) with a trivial form factor, k = 1 (dashed yellow) and 
k = 3 (dashed orange) are shown. For all lines, we use, Vq = 220 km/sec, fesc = 500 km/sec, 
and take the Earth velocity on June 2nd. 



2.1 Dark Matter velocity distributions 

The DM interaction rate with nuclei depend on the velocity distribution /^ of DM as in 
eq. (2.1). Since this distribution has not been measured, it adds its own source of uncertainty 
to the DM interpretation of direct detection experiments. This topic has been extensively 
studied in the literature, where several possibilities such as streams or velocity substructure 
were considered [23, 27, 28, 29]. While a systematic study of astrophysics uncertainties and 
its prospects for ameliorating the experimental tension is beyond the scope of this paper, we 
consider few motivated distributions to demonstrate the possible variation in the fits. We 
stress that it is conceivable for significantly different conclusions to be drawn with the use of 
less conventional or more exotic possibilities, as was demonstrated prior to the new CoGeNT 
result in, e.g. [27]. 

The velocity distribution in our local frame, f®{v,t) defined in Eq. (2.1), is conveniently 
expressed in terms of the velocity distribution in the galactic frame, /(f), through 



f®iv,t) = f{v + V(Q{t);vo,Ve 



(2.6) 



Here % is the relative motion of the Earth with respect to the galactic frame (see e.g. [28] 
for further details). Vq is the root mean square velocity typically taken to be in the range. 



220 < f < 270 km/s, while fesc is the escape velocity, in the range, 450 < fcsc < 650 km/s [22]. 
Below, we consider vq = (220, 270) km/s and fesc = (500, 600) km/s. 

The DM velocity distribution in the galactic frame is often assumed to be a Maxwell- 
Boltzmann (MB) sharply cut off by a finite escape velocity, 

/mb(w; ^^0, ^esc) = ^rre~'''^''°0(^esc " v) , (2.7) 

i\e 

with Ne = (erf(2;)— 2zexp(— 2;^)7r~-'^/^)7r^/^t>Q and z = v^sc/vq. Above, we have explicitly denoted 
the dependence of the velocity distribution on vq and fesc- The above MB distribution does not 
seem to capture the results of A^-body simulations [24, 25, 37, 38]. An improved ansatz for an 
isotropic velocity distribution is given by 
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for 1.5 < k < 3.5 [29]. The MB distribution is reobtained in the limit A; — ?■ 0. These velocity 
distributions are plotted in Fig. la. 

The above distributions enter the scattering rate, Eq. (2.1), through the function,^ 

""" "^-^^^ '' dhUv,t)(^] ^FUq,v). (2.9) 



^^0 Vo Vo ) J\v\>v^,,, V^O. 

Note that the dependence on V(q and fesc is implicit through the definition of /® [see Eq. (2.6)]. 
In Fig. lb we plot ?7„ for F'^y^i^q.v) = 1,v'^/vq,v'^/vq with the Maxwell-Boltzmann distribution, 
as well as the trivial form factor for the fk velocity distribution defined in Eq. (2.8). As can 
be seen, higher powers of velocity in the form factor imply a larger overall normalization and a 
larger sensitivity to higher minimal velocity with a sharper falloff. We return to these features 
in r]n in Section 4.7, where we study velocity-dependent form factors. 

3 Experiments 

In this section we briefly summarize the data used to derive the allowed region for spin- 
independent DM scattering. For each experiment we stress the various sources of uncertainties 
and our approach for dealing with them. In most cases, we consider several possibilities, allow- 
ing for a conservative view of the tension between the positive and null experimental results. 

3.1 DAMA 

The DAMA experiment employs a Nal(Tl) target and observes an 8.9o" evidence for an annual 
modulation in its energy spectrum [1]. The modulation is present in the 2-6 keVee energy 
range and the time dependence of the rate is consistent with the hypothesis of DM scattering 



^We use the notations of [28], however our function rjn is chosen to be dimcnsionless, as opposed to the 
corresponding function defined in [28]. 
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2-2.5 


0.016 


0.004 


2.5-3 


0.026 


0.005 


3-3.5 


0.022 


0.005 


3.5-4 


0.008 


0.005 


4-4.5 


0.011 


0.004 


4.5-5 


0.005 


0.004 


5-5.5 


0.009 


0.003 


5.5-6 


0.004 


0.003 


6-14 


0.000 


0.000 



Table 1: The data used in the fit to the DAM A modulated amplitude. 

on nuclei. The DAM A experiment is able to measure only the fraction of energy that recoil 
nuclei deposit as scintillation. This fraction of the total recoil energy, known as the quenching 
factor, is taken to be qi = 0.09 for Iodine. 

Different groups report varying values for the quenching factor on Sodium. The DAMA 
collaboration reports gNa = 0.30 ± 0.01 averaging over recoil energies ranging from 6.5 to 97 
keV. Ref. [39] finds 0.33 ±0.15 between 4 to 11 keV while ref. [40] finds 0.252 ±0.065 around 10 
keV. The relevance of gNa resides in the fact that larger values ameliorate the apparent tension 
between DAMA and the null experiments. Indeed larger gNa implies lower recoil energies at 
DAMA and consequently favors smaller WIMP masses. Following many previous works, we 
adopt a conservative estimation of the uncertainties, assuming gNa = 0.3 ± 0.1. 

An additional source of uncertainty arises due to the crystalline nature of the target material. 
It is possible that some (experimentally unknown) fraction of the ion scatterings occur parallel 
to a symmetry axis (channeled events), depositing the entire energy in scintillation. For such 
events the quenching factor is effectively 1 and hence a sizable fraction of channelled events may 
significantly alter the direct detection predictions [31]. Recent theoretical results, suggest that 
the channeling fraction in Nal is negligible [33]. In order to be over-conservative, in Section 4.9 
we also perform fits in which we allow an energy independent channeling fraction. 

The spin independent fit to the DAMA modulated rate allows two qualitatively different 
best-fit regions: one around M^m ~ 80 GeV with a ~ 10~^^ cm^ due to scattering on iodine 
{A = 127, Z = 53) and one around M^m ~ 10 GeV with a ^ 10~^° cm^ due to scattering on 
sodium [A = 23, Z = 11). The former region is firmly excluded by many other experiments, 
most notably XENONlO/100. The latter region, while still disfavored by other null searches, is 
not as badly excluded due to many experimental uncertainties and due to the general difficulty 
of direct detection searches to deal with low recoil energy scatterings. In the following we will 
focus on this low mass region, being also the only possibility to reconcile DAMA with the 
results of CoGeNT. 
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^3As 


12.7 


0.33 


1.414 


0.078 


80 


68Ge 
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0.01 


1.298 


0.077 
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52.8 


0.05 


1.194 


0.076 


271 


65Zn 


211 


0.02 


1.096 


0.076 
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56Ni 
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0.23 


0.926 
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0.45 


0.846 
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0.846 


0.074 
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44.9 


0.12 


0.769 


0.074 


996 


5^Mn 


21.1 


0.09 


0.695 


0.074 


312 
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2.94 


0.15 


0.628 


0.073 


28 


49y 


14.9 


0.12 


0.564 


0.073 


330 



Table 2: The data used to extract the L-shell backgrounds. N is the total number of decays 
expected from a given isotope from the beginning of the CoGeNT data taking to the end of 
time, with 6N/N being its relative error. -Epeak is the central value for the corresponding binding 
energy while aE ^^^^ is the energy resolution. T1/2 is the half-life of the relevant isotope. 

To fit the modulated signal at DAMA we build a simple x^ using the content of the first 8 
bins in Fig. 9 of [1]. We use a single bin for the modulated rate in the energy range of 6 to 
14 keVee. The data we use are shown in Table 1. DAMA collaboration also reports (see fig. 1 
of [41]) a measurement of its non-modulated rate. This can be used to define an upper bound 
on the signal, following for instance the procedure of [42]. Unless stated otherwise, we will not 
use this bound in the fits. 

To compare the data with the theoretical hypothesis we introduce a finite detector resolution, 
parameterizing the energy smearing through a gaussian with energy dependent width given by. 



C"DAMA 



=0.448W- + 9.1x10^3- . (3.1) 

keVee/ V keVee keVee 



3.2 CoGeNT 



The CoGeNT experiment uses a Germanium {A ^ 76, Z = 32) detector and takes data in 
the Soudan Underground Laboratory (SUL). We employ the latest data release [3], which are 
the result of 442 live days of data-taking from January 4 2010 to March 6 2011 on a 0.33 kg 
Germanium target. The new data confirm the presence of an exponential distribution of events 
between 0.5 and 1.5 keV which is not accounted for by any known background. A time analysis 
of the same data also shows evidence of an annual modulation in the 0.5-3.0 keV range, which 
could be interpreted as evidence for DM interacting with the detector. 

To extract the rate and the modulation to be used in the fit, we use the time-stamped 
raw data obtained from the CoGeNT collaboration [43]. The relevant region between 0.4 and 
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Figure 2: Left: distribution of the events observed during 44^ days of live-data taking (black 
dashed). The black line is the result of efficiency unfolding. The blue contour shows the contribu- 
tion from L-shell EC lines calculated from the data in tab. 2. Right: CoGeNT non-modulated 
rate after the subtraction of L-shell lines. The blue line is the result of an exponential+ constant 
fit to the rate. In both plots we use a uniform (0.1 keVee) binning. 



3.5 keV is contaminated by events due to the electron capture (EC) decay of cosmogenicaUy 
activated elements in the detector. In a given time window these events are expected to show 
up as peaks centered around the L-shell binding energy of the daughter of the decaying nucleus. 
The width of the peaks is related to the detector resolution and it's amplitude to the number 
of active isotopes in the detector. The amplitude of the peak is thus expected to decay in time 
with the half-life of the relative isotope. In Table 2 we show the parameters needed to describe 
these background events. It is clear that due to the time dependence of such backgrounds, their 
subtraction is of primary importance before attempting any time-analysis of the signal. 

In Fig. 2 we show the spectrum of events observed at CoGeNT as well as the rate after 
efficiency unfolding. As shown, the residual events after subtracting the L-shell peaks, are well 
fitted by an exponential plus constant function Ae~^^ + C, with E in keVee. We find. 



A = 100.4, B = 3.4, C = 2.4, xVd-o.f . = 57/(50 - 3) . 



(3.2) 



For the time-analysis we use data up to 6 keVee binned into 5 energy bins: 0.5-1, 1-2, 2-3, 
3-4.5, 4.5-6 keVee. Each bin is further divided into 15 time bins of 30 days each. For every 
energy-time bin we subtract the estimated number of events from L-shell EC and apply a 
correction to account for the offiine time. Since in the first bin the efficiency is not constant, 
we choose not to unfold it to avoid loosing more statistics. A first test to the DM hypothesis 
is to check whether the signal modulates with a yearly period peaking on 2'^'^ of June. We thus 

fit the signal to the function 

27r(t - (f)) 



B + S cos 



T 



(3.3) 



and check the consistency with the hypothesis T = 365 days and = 152 days with the use of 
a x^ analysis, marginalizing over B and S. The results of the fit are shown in Fig. 3. The data 
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Figure 3: Best fit for the period T and phase (f) of the CoG^l^T signal in the low energy bins. 
The continuous lines are the la contours. For the large 0.5-3 keVee bin, we also show the 2a 
one (dashed). The black cross indicates the prediction of the DM hypothesis. 

show consistency with the DM hypothesis at the 2a level. The signal deviates from the null 
(no-modulation) hypothesis, at the 2a level at most for each bin separately, while it exceeds 
this value (however not reaching 3-sigma) for the combined 0.5-3 keVee bin. The deviation 
from the null hypothesis is demonstrated further in Fig. 5, where the modulation spectrum is 
shown. 

We thus fix T to the value predicted by the DM hypothesis, 365 days, and proceed to 
determine the best fit to the other parameters B, S and 0. The results are summarized in 
Table 3 and shown in Fig. 4. The dashed curves in Fig. 4 are obtained by fixing both T = 365 
days and (j) = 152 days. To extract the amplitude of the modulation S to be used in what 
follows, we marginalize over B to obtain the spectrum shown in Fig. 5. In the two bins 0.5 — 3 
and 3 — 6 keVee we get (in cpd^^kg^^keVee"^) 

0.5-3keVee: 0.43±0.18, 3 - 6keVee : 0.02 ± 0.11 (3.4) 

Finally, we need to specify the quenching factor of Germanium in order to translate keVee 
energy into keV nuclear recoil energies. We use a Lindhard k = 2 parametrization 

^■A^r) ■ (3-5) 



keVee V keV 



3.3 SIMPLE 

The SIMPLE experiment is a superheated liquid C2CI F5 droplet detector. Working in a manner 
similar to bubble chambers it looks for bubble nucleations induced by WIMPs. In particular, 
due to the presence of light ions such as fluorine {A = 19, Z = 9), it is sensitive to light 
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Figure 4: Time- analysis of the data in the bins of interest. In all the fits we keep T = 365 days 
fixed. The blue full lines are obtained floating B, S and cf). The blue dashed lines are the results 
of the fit fixing also cf) = 152 days. 



E (keVee) 



E 



4> 



x^ 



0.5-1.0 47.3(47.6) 4.6(3.4) 

1.0-2.0 26.1(26.3) 6.4(4.3) 

2.0-3.0 22.6(22.8) 3.5(2.1) 

3.0-4.5 34.2(34.2) 0.2(0.1) 

4.5-6.0 46.9(47.0) 1.9(0.3) 



114(152) 7.9(9.1) 

108(152) 11.5(13.7) 

104(152) 8.0(10.1) 

211(152) 15.7(15.7) 

77(152) 21.2(21.7) 



0.5-3.0 97.2(97.7) 15.0(10.7) 112(152) 6.3(11.0) 
3.0-6.0 82.5(82.7) 2.7(0.7) 81(152) 14.3(14.9) 

Table 3: Summary of the time-analysis of the CoGeNT data. Everywhere T = 365 is fixed. 
The numbers in parenthesis are obtained fixing (p = 152 days. The number of d.o.f. are 12 
(13), in each energy bin. 

mass WIMPs even with a relatively high threshold of 8 keV. Although it is not possible for 
this kind of experimental setup to measure events' energies, the threshold can be set precisely 
as the bubble nucleation depends on the temperature and pressure of the superheated liquid. 
Neutron induced recoils has confirmed the minimum threshold energy at 8 keV with a precision 
of O.lkeV. 

The SIMPLE collaboration has recently published new results [44] including new Stage 2 
data which improves and merges the analysis of older (Stage 1) data. As the merging can 
potentially introduce additional systematic uncertainties, we use only the Stage 2 data, with 
zero unidentified events and an exposure of 6.71 kg days. We set an exclusion limit using a 
Poissonian likelihood 
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-A''dm 



(3.6) 



Such a likelihood is considered to be a conservative choice that produces bounds compatible 
with the official one. 
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Figure 5: Spectrum of the modulation amplitude in the CoGeNT data assuming T = 365 days 
and (f) = 152 days. In red we show the result of fitting just two bins: 0.5-3 and 3-6 keVee. No 
efficiency correction is applied. 

3.4 XENONIOO 

XenonIOO is a two-phase Xenon {A ^ 131, Z = 54) experiment which pubhshed results 
obtained from approximately 100 live days of data acquisition in a fiducial volume of 48 kg. The 
consistency of the outcome with the background hypothesis allows to place strong constraint 
on the interaction of a WIMP with Xenon nuclei. 

As opposed to DAMA and CoGeNT, XenonIOO has a signal-to-background discrimination 
ability. This is achieved by comparing the primary scintillation signal (SI) to the ionization 
yield (S2), the relative magnitude of the latter being bigger for electronic recoils. The capability 
of XenonIOO to detect low mass WIMP scatterings crucially depends on the response function 
£eff which, through the relation 



SliEn) = 3.6 FExErX C,s, 



(3.7) 



gives the number of photoelectrons (PE) in the SI signal as a function of the recoil energy. 
Measurements of Ces extend down to 3 keV [45] while extrapolation have to be used for lower 
values. It was argued in [34], that the uncertainty on C^s significantly influences the ability of 
XenonIOO to constrain light DM. Below, we use the £es contours adopted by the XenonIOO 
collaboration and shown in Fig. 1 of [6]. To have a good signal-to-background discrimination, 
the lower XenonIOO threshold is fixed at 4 photoelectrons which corresponds to roughly 8 keV 
(depending on the precise choice of Ces)- The relevance of the extrapolation of Ces to lower 
energies has to do with the statistical nature of the scintillation process. Assuming the photo- 
electron generation process to be poissonian in nature, with a mean dictated by Eq. (3.7), recoil 
occurring below threshold will have a non vanishing probability to generate an SI signal above 
threshold. This tail is crucial to the constraining power of XenonIOO for low mass WIMPs^. 



^Thc XenonIOO analysis considers an S2 threshold of 300 photoelectrons. Values of (SI) below threshold 
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Figure 6: Xenon100/10 bound compared with the 3a favored region for DAMA (magenta) 
and CoGeNT (red). The blue band is the XenonIOO exclusion obtained by varying £eff in 
its 1 o range. The blue dashed line shows the result of the conservative choice Ccs = below 
4keV. The light-green band represents how the 3a XenonIO exclusion changes with respect to 
Qy variations: on the left we use the choice of the XenonIO collaboration and on the right we 
use the conservative estimate adopted in [34]- Two choices for the sodium quenching factor q^a 
are adopted for the DAMA fits. 



We have checked and found that assuming a different statistical behaviour of the Sl-generation 
process (binomial for instance, see [34]) does not alter our conclusions. We show the XenonIOO 
bound and its uncertainty in Fig. 6. The blue band is the 3o" exclusion obtained by varying £eff 
in its la range [6]. The blue dashed line comes from the conservative choice of having £cff = 
below 4 keV. 

To determine the XenonIOO bounds we follow [47] doing an event-by-event fit to the three 
observed events at energies 6, 19 and 22 keV assuming a uniform background in the region 
4PE < SI < 30 PE normalized to the total number of expected events, 1.8. We include finite 
resolution effects through a gaussian smearing of the SI signal with a width given by O.SvSl. 
We include an SI peak finding efficiency in the same way it is done in [46]. We find this to be 
a small correction. 



3.5 XENONIO 

The high XenonIOO (XenonIO) threshold, roughly 8 keV (5 keV), is related to the fact 
that recoils of too small energy are not efficiently converted into a primary scintillation signal 
SI. It is through comparison of this signal to the secondary scintillation (S2 signal) that the 
background from electronic recoil is subtracted. 



may generate, through a statistical fluctuation, an SI signal above threshold but may fail to pass the S2 cut. 
Following [46] we ignore recoils giving (SI) < 1 PE. 
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A low threshold analysis (1 keV) of Xenon 10 data is available [48] by discarding the SI 
signal altogether and looking just at the S2 signal. Backgrounds are thus allowed to pollute 
the data but the gain in sensitivity at lower masses is substantial. 

Since at low recoil energies the scintillation photons do not give a measurable SI signal, one 
is forced to calibrate the energy scale using just the S2 signal. This calibration is encoded in 
the so called Qy parameter defined as 

Q.(Er) - ^, (3.8) 

where N^ is the number of measured photoelectrons for a recoil of a given energy E^i. The role 
of Qy is similar to the one of C^s- Fixing a lower threshold for the magnitude of the S2 signal 
will, depending on Qy, define a lower threshold on the nuclear recoil energies probed by the 
experiment: the smaller the Qy the higher the effective threshold. Since the value of Qy is not 
measured below 4 keV [49, 50] extrapolations must be used. 

XenonIO takes its S2 threshold at 5 photoelectrons, corresponding to roughly 1.4 keV with 
their choice of Qy. According to [34] this choice is far too generous. [34] thus proposes a smaller 
Qy (see details in the original reference) corresponding to an energy threshold of roughly 4 keV. 

To calculate the XenonIO bounds we use the observed events in Fig. 3 of [48], corresponding 
to 15 kg days of effective exposure, passing all the 5 cuts detailed in Tab. 1 therein. In a fashion 
which resemble the Pmax method [51], we determine the two consecutive events between which 
the total expected signal from DM is maximized. We use poissonian statistic to define the 
likelihood, 

-Cxeio = e-^'™ , (3.9) 

where S2 is the number of photoelectrons (including acceptances and poissonian statistical 
fluctuation which affects the S2 resolution) expected in the above interval under the DM hy- 
pothesis. We thus use Xxeio — ~2 In £xeio to set bounds. In the rest of the paper our choice 
for Qy is the same as the one adopted by the XenonIO collaboration. A comparison with the 
more conservative choice of [34] is shown in Fig. 6. 

3.6 CDMS-Si 

The cryogenic CDMS experiment ( performed at the SUL like CDMS) operates Silicon {A ^ 28, 
Z = 14) and Germanium solid-sate detectors. Like XenonIOO, CDMS has the ability to 
discriminate between nuclear recoils and electronic backgrounds, measuring both ionization 
and phonon signals. We use the results from the unofficial analysis [52, 53] of 6 Si detectors 
corresponding to a raw exposure of 53.5 kgd. We use the efficiency reported in [52, 53] which 
drop to zero below 7.8 keV. No events are observed with an expectation of 1.1 events from 
surface backgrounds. We use the poissonian likelihood, 

£gi = e"^™ , (3.10) 

where Ndm is the expected number of events under the DM hypothesis to extract the bounds. 
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As noted in [54], CDMS-Si data indicate a behavior for the Sihcon quenching factor which 
is not in accord with the prediction of Lindhard theory (see for instance Fig. 3.20 of [52]). This 
can be amehorated by correcting the energy scale by (9(20%), which also goes in the direction 
of weakening the CDMS-Si bounds but we do not include this correction in the fits. 

3.7 CDMS-Ge 

For the analysis of CDMS Germanium we employ the data obtained in the recent low-energy 
analysis [55], with an energy threshold of 2 keV. Since the electronic and the nuclear recoil 
bands merge at these low energies, the sensitivity to low WIMP masses offered by this low- 
threshold analysis comes at the price of accepting a large amount of background. Though the 
CDMS collaboration provides possible explanations for the background and claims that it can 
explain all the events in the signal region, they neglect it as being dependent on extrapolation 
and on too many uncertainties. This approach has been recently criticized in various ways [56] . 
In our analysis we follow the approach held by the CDMS collaboration not including back- 
grounds and treating, conservatively, all observed events as DM induced recoils. We depart, 
for simplicity, from the CDMS group retaining only the data coming from the best performing 
detector (T1Z5), that observed 36 events between 2 and 20 keV. By considering all the events 
as signal we construct a x^ by fitting the theoretical total rate A^dm, so that 

2 (A'dm - 36) /Qii\ 

XGe = 7^ e(A^DM-36). (3.11) 

Our bound agrees with the one reported by the CDMS collaboration. Due to the large number 
of observed events (and in the absence of background subtraction), the x^ function in Eq. 
(3.11) is rapidly growing thus giving confidence levels which are much closer to each other in 
comparison to the other experiments. This causes a deterioration of the global fit once the 
CDMS-Ge bound is included. 

4 Results 

In this section we present the best fit results for the experiments discussed above. Our approach 
here is to study the influence of different forms of cross-sections and velocity distributions, as 
discussed in the introduction, on the fits to the data. This approach differs somewhat from the 
usual effective theory one, where the bounds on different operators are studied. The virtue of 
the current method is that it allows one to identify the physical necessities in order to minimize 
the experimental tension. Realistic scenarios of course, often require taking linear combinations 
in either approaches. 

4.1 Format of the figures 

In our figures, all continuous (dotted) curves correspond to 95% (99.7%) C.L. for two degrees 
of freedom, i.e. 2a (3o") corresponding to Ax^ = 6 (Ax^ = 11.6). In the figures we show the 
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following: 

• In purple, the region favored by the modulation observed by DAMA at fixed value of gNa- 
The parameter fixed by the fit is essentially Mdm x Q'Na, such that considering a higher 
gNa value shifts the favored region to lower DM masses. The text 'DAMA' lies over the 
DAMA best fit point. 

• In yellow, the wide region favored by the modulation observed by CoGeNT. We further 
impose that the DM rate alone does not exceed the rate observed by CoGeNT in any 
point. In view of the poor statistical significance, we plot here the 68% CL. contour 
(dashed-dotted) as well as the 95% CL. contour. 

• In red, the favored region for the rate observed by CoGeNT (data in the right panel of 
fig. 2) assuming, in addition to the DM signal, an energy-independent constant term with 
normalization fixed as in Eq. (3.2). This assumes, arbitrarily, that the L-shell decays of 
activated isotopes account for the majority of low energy background events. The region 
favored by the CoGeNT rate up to this caveat is shown in red and marked as 'CR'. 

• A green curve for the bound obtained with the CDMS Silicon result, denoted as 'Si'. 

• A red curve for the bound by CDMS with Germanium, denoted as 'Ge'. 

• A blue curve for the bound by XenonIOO, denoted as 'Xeioo'- 

• A purple curve for the bound by XenonIO, denoted as 'Xeio'. 

• A dark yellow curve for the bound from SIMPLE, denoted as 'C CI F'. 

• Finally, we plot a green dot indicating the global best fit taking into account all signals 
and bounds. 

We define the global ~x^ as x^ = ^j X?) summing over all the relevant experiments, and report 
the value of x^ for the best fit point. Such x^ can be used as statistical indicator to compare 
different fits and to evaluate the overall quality of the fit (although other more sensitive statisti- 
cal indicators exist). For experiments that give bounds we fixed x? = when no DM is present, 
such that we expect that a good fit should correspond to y^ ~ nobs — ""-par, where nobs is the 
number of observed data-points and npar the number of free parameters. Since the CoGeNT 
rate could be contaminated by unknown backgrounds at low energy, we perform two global fits: 

i) Fitting all data including the CoGeNT rate. A good fit should have xf) ~ 43, as we fit 
35 data points in the CoGeNT rate, 9 in the DAMA modulation, 2 in the CoGeNT 
modulation with a number of free parameters going from 2 to 4; 

ii) Fitting all data (the CoGeNT modulation and all other experiments) but dropping the 



CoGeNT rate. A good fit should have xfi) ~ 8. 

In the pictures we report the value of both y^ evaluated at their best-fits, using the format 
"x^ = Xn^xfi)"- O^^ Results are presented in the next sub-sections. 
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Figure 7: Upper left: Standard fit. DAMA and CoGeNT do not overlap, and are excluded 
by many experiments. In the other plots we vary the DM velocity distribution, finding minor 
changes. On the top right: we use the smooth distribution of Eq. (2.8) with k = 3. Bottom 
left A higher Vq = 270 km/s and lower v^sc = 500kni/s are assumed. Bottom right: A higher 
I'esc = 600 km/s is taken. In all plots fp/fn = 1 d'^d q^a = 0.3. See Section 4-1 for the color 
coding. 

4.2 Standard fit 

In Fig. 7a we show the "standard" fit, in terms of elastic spin- independent DM, using the cross- 
section in Eq. (2.2), and assuming no form factor, Fdm = 1 and fn = fp = 1- We see that (i) 
DAMA and CoGeNT do not overlap and (ii) they are excluded or strongly disfavored by many 
experiments. As a result the global best fit (green dot) has a very high x^, and corresponds to 
roughly no effect in DAMA. 
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Figure 8: Top left: Standard fit assuming a higher quenching factor for Sodium, gNa = 0.4. 
Top right: Global fit assuming is o spin-violation, dictated by the fp/fn parameter. The sit- 
uation significantly improves, but the best fit (marked with a green dot) remains very poor. 
Color coding is described in Section 4-1- Bottom: DM predictions for the best fit point, al- 
lowing a floating fp/fn- The signal is plotted against the DAMA and CoGeNT modulated and 
unmodulated data. 

4.3 Astrophysical uncertainties 

We explore the sensitivity to modifications in the velocity distributions (under the assumption 
of isotropy) by first considering the smoothed cuts discussed in section 2.1, as controlled by 
the parameter k: bigger k implies a smoother distribution while the sharply cut, Maxwell 
Boltzmann distribution is obtained in the fc — )■ limit. These velocity distributions are shown 
in Fig. la. As can be seen in Fig. 7b the fits to the experimental data assuming a DM velocity 
distribution with a smooth k = 3 cut are quite similar to the ones for k -^ (Fig. 7a). In 
view of the very minor difference from now on we stick to sharp cuts and we do not show the 
intermediate case k = 1. 

The improvement in changing vq is also small. The slope of DAMA (CoGeNT) spectrum 
fixes the value of fi'^VQ/rriN, hence for low-mass WIMPs raising vq favors smaller DM masses. 
The effects are shown in Fig. 7c. Similarly, changing the maximal DM velocity fgsc does not 
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Figure 9: Ranges of transferred momentum q (left) and minimal DM velocities, fmin; (right) 
probed by relevant experiments. Solid bands: CoGeNT 0.5-3.0 keVee band (yellow) and DAMA 
2-6 keVee bands (purple) for both unchanneled Sodium (upper panel) and channeled Iodine 
(lower panel). Dashed lines: XENON 100 3 keV threshold (blue), CDMS-Si 7 keV threshold 
(green), CDMS-Ge 2 keV threshold (red) and SIMPLE 8 keV threshold (dark yellow). The 
horizontal lines show the maximum possible DM velocity in the Earth frame for fesc = 500 
km/s (thick) and fesc = 600 km/s (thin). 

help, see Fig. 7d. 

4.4 Quenching factor uncertainties in DAMA 

Assuming a quenching factor for Na higher than what is claimed by the DAMA collaboration, 
allows to shift the DAMA best fit region to lower Mdm (the combination Mdm x Q'Na being 
essentially fixed), improving the fit. In Fig. 8a we show the case gNa = 0.4, that we will adopt 
from now on. Even so, the best fit is so bad that we do not show its comparison with data. 

4.5 Isospin violating couplings 

It is possible that DM does not couple equally to protons and neutrons. Such isospin- violating 
DM scattering has been studied against the CoGeNT rate data previously in [15, 16, 17]. 
In Fig. 8b (and in the rest of the fits) we relax the assumption of fp = fn, and allow for 
different spin-independent cross sections on protons and neutrons. As can be seen, fioating the 
additional parameter, fp/fn, allows to 

• Improve the agreement between CoGeNT and DAMA. 

• Ameliorate the tension of the positive results with one experimental bound, through the 
tuning of fpZ + fn{A — Z) ^ 0. In our fits, we take into account the different isotopes in 
the various experiments, thereby allowing for only a partial cancelation. 

When performing a global fit, we find the best-fit shown in Fig. 8b, with a significant improve- 
ment over the fp = fn case shown in Fig. 8a. Still, the global fit remains very poor, despite the 
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Figure 10: Fit with free fp/fn and with 
momentum- dependent form factor. In the 
left upper (lower) row the cross section 
grows (decreases) with transferred momen- 
tum q^. The top right plot shows the fit for 
a q^ form factor. Color coding is given in 
section 4-1- 



assumption of a relatively large Sodium quenching factor, gNa = 0.4. Indeed, the incompatibil- 
ity of the tentative DM signals is with four experiments, performed M^ith three different nuclei 
(Xe, Si and Ge). Consequently, this extra parameter allows to avoid the strongest constraint 
(coming from Xe experiments), but is insufficient to relieve the tension completely. 

4.6 Momentum-dependent elastic scattering 

Momentum dependent scattering arise in several instances. One notable case is when the DM- 
nucleon interaction is mediated by a pseudo-scalar. It is possible to systematically study the 
momentum dependent effects by utilizing non-relativistic effective theory. Indeed, the typical 
energy transfer in DM-nucleus collisions relevant for direct-detection is much below the nuclear 
binding energy. Additionally, the DM velocity, v ~ 10~^, is non-relativistic. Consequently, it 
is possible to describe the scatterings via an effective theory assuming a rotationally invariant 
potential [57]. Restricting to spin-independent interactions, one finds to leading order in the 



20 



DM velocity, f , and the moinentuni transfer, g, only four scalar operators from which the 
effective potential may be constructed. Written in momentum space they are, 

^1 oc 1, a (X 1 , 

M oc Sdm ■ g, cr oc g^ 

-4-3 oc Sdm ■ "y , cr oc t;^ , 

^4 oc s ■ g X t7, cr oc g^f ^ , 

where sdm the spin of the DM. In a Lorentz invariant theory, specific linear combinations of 
the ^s appear in the full DMA^ -^ DMA^ amplitude. For further discussion, see [57]. 

From a phenomenological point of view we can study the effect of v- and g-dependent form 
factor on the compatibility between DAMA, CoGeNT and the null-experiments. We consider 
the following DM form factors that enter in Eq. (2.2): 

where the momentum transfer is normalized with respect to the reference gref = 100 MeV. The 
last form factor can be generated in models where DM-nucleon scattering is mediated by a new 
particle with mass rti <^ gref. The 1/g^ form factor can arise in this model with a coupling of 
type A2- 

A momentum-dependent form factor changes the relative rate in different experiments ac- 
cording to the range of g^ they probe, as summarized in Fig. 9a. We show bands corresponding 
to the signal regions for DAMA (2-6 keVee) and for CoGeNT (0.5-3 keVee), as well as lower 
g^ value corresponding to the energy threshold for all other experiments. For SIMPLE we show 
the g^ corresponding to fluorine which is target element with the largest rate thus being the 
main contributor to bounds. 

g^ and g^ form factors help in this respect since they deplete the spectrum at low recoil 
energies and hence require smaller masses for a good fit. The Xenon 100 bounds are expected 
to get stronger relative to the other experiment due to the higher momentum transfer. It must 
be noted, however, that in the low mass region, the sensitivity of XenonIOO is saturated by 
its threshold (and consequently the bound becomes vertical at small masses). Another feature 
of g" in > 0) form factors is that they enhance the CoGeNT signal with respect to the DAMA 
one, since the former probes a higher g region. This implies that the ratio of the best-fit cross 
sections for CoGeNT and DAMA will be smaller with respect to the standard fit, making the 
separation of the two regions more severe. This problem can be ameliorated by an appropriate 
choice of /„ 7^ fp. Fig. 10a shows the global fit. We find that g^ form factor gives mildly better 
fits than the standard q^ case, however the overall x^ remains poor. Similar results are found 
for a g^ form factor. 

An analog, but inverted, discussion apply for 1/g^ or 1/g^ form factors. These form factors 
move both the DAMA and CoGeNT regions to higher masses and make the two best-fit points 
closer (assuming /„ = fp). The DAMA non-modulated signal, which probes recoil energies 
below the ones where modulation is observed, excludes the g~^ behavior. Fig. 10b shows the 
results for 1/g^. We do not find that the form factor improves the global fit. 
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Figure 11: Fit with free fp/fn cind with a velocity dependent form factor. Left: f^ form factor. 
Right: f^ form factor. An overall normalization effect is visible, as discussed in Section 2.1. 
For further discussion see Section J^.l. Color coding is given in section J^.l. 

4.7 Velocity-dependent elastic scattering 

The amphtudes in Eq. (4.1) also induce velocity-dependent form factors. Here we take again 
the phenomenological approach and study the effect of the form-factors, 



p2 



2 / 2 



P2 



4 / 4 

V jv^ 



(4.3) 



where we normalize with respect to a mean velocity Vq. 

The effect of the overall normalization discussed in Section 2.1 and shown in Fig. 1 is trivial, 
as can be seen in the fits displayed in Fig. 11. To understand the effect of the sensitivity to 
the minimal velocity, we show the minimal DM velocity probed by different experiments as 
a function of the DM mass in Fig. 9b. We note that the minimal velocity does not depend 
on astrophysics (i.e. velocity distributions) nor on the scattering cross sections (see also [26]). 
The bands for CoGeNT and DAMA corresponding to the observed signal energy range are 
shown. For null-experiments we show the fmm corresponding to the respective energy threshold, 
where we have taken 3 keV as a representative value for XenonIOO. For completeness, the 
horizontal lines show the maximum velocity that a DM can have with respect to the Earth, 
given fesc = 500 or 600 km/s. We see that neither CoGeNT not DAMA probe values of v 
inaccessible to other experiments. 

Nonetheless, some improvement in the global fit occurs when the typical minimal velocity 
is around the tail of the ?7„ function. In that region, a small change in sensitivity to f mm results 
in a large change in the scattering rate. Indeed our global fits shown in Fig. 11 demonstrate a 
small improvement in the fits. 
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Figure 12: Fit with inelastic dark matter and free fp/fn o-nd DM splitting, 5. Left: An example 
of endothermic DM (6 > 0), which improve the global fit. Right: exothermic DM (6 < 0) 
improves the partial fit. Color coding is described in Section 4-1- Bottom: DM predictions at 
the best fit point in the right panel (exothermic DM), plotted against the DAMA and CoGeNT 
modulated and unmodulated data. 

4.8 Inelastic light Dark Matter scattering 

It is possible that DM scatters inelastically with the nucleus. For inelastic scattering to take 
place, two semi-degenerate DM states are assumed, with mass splitting 6 = M^j^ — Mdm- Up- 
scattering of the lighter DM state requires it to have enough energy, thereby suppressing the 
rate for small values of the recoil energy. Since up-scattered iDM kinematically favors heavy 
targets, it was originally able to ameliorate the tension between the DAMA modulation and 
the null CDMS resuh [9], requiring Mdm ~ 100 GeV and 6 ~ 100 keV. 

This original iDM scenario seems now excluded by the recent XenonIOO results [58, 59]. 
Here we study a different iDM regime with significantly smaller splitting and lighter DM (for 
previous study see [15]). In this window the tension with the null experiments is ameliorated 
mostly due to the small DM mass, while the DM scattering rate is falling above a few keV, 
much like in the elastic case. The effect of 6, is to modify the minimal velocity needed for 
scattering to occur, in accord with Eq. (2.4). 
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For the masses and recoil energies of interest, the sphtting 6 is required to be smaller than 
~ 15 keV in order to comply with the CoGeNT data. Consequently, the inelasticity for 
up-scattering is only relevant for the light Sodium and Silicon targets in DAMA and CDMS 
respectively. We find that such scenario does improve the global fit only slightly, and one 
example for 6 ^ 11 keV is shown in Fig. 12a 

It is possible, however, that the heavier DM state is cosmologically long lived, in which case 
it occupies a sizable fraction of the DM density. In fact, this possibility occurs quite naturally, 
as noted in [60, 61]. In such a case, DM can also down-scatter with the nucleus, producing an 
exothermic reaction. This possibility, dubbed exothermic DM (exoDM), was studied in [61]. 
Technically, the rate for exoDM is given by the same expression as for the up-scattered iDM 
case, with now 6 < 0. But the behavior of exoDM is significantly different. The minimal 
velocity, Eq. (2.4), is minimized for Eji ^ \6\^/mN and hence lighter targets are more sensitive 
to exothermic DM scatterings. No net modulation in the total number of signal events is 
expected. However, since the spread of the spectrum around its maximum depends on the DM 
kinetic energy, the recoil energy spectrum will actually modulate annually. The modulated 
spectrum observed by DAMA can thus be reproduced [61]. 

In Fig. 12b we show the best fit for exoDM. In order to isolate the effect of down-scattering, 
here we only consider down-scattered events, even though up-scattering is expected to be sig- 
nificant. For the fit, we fioat the DM mass, the cross section, fp/fn and the DM mass splitting, 
finding that a (5 ~ —10 keV improves the partial fit with respect to the analogous fit with 6 = 0, 
shown in Fig. 8b. 

4.9 Channeling 

We now allow for some amount of channeling in the Nal(Tl) crystals of DAMA and/or in 
CoGeNT. In this section we do not introduce any inelasticity: 6 = 0. Fig. 13a shows how 
a relatively small fraction of channeled events is enough to shift the best fit regions for these 
experiments to a different global minimum of the x^- Next, increasing the channeling fraction 
further, mildly shifts the best-fit regions down to lower values of the cross section, but no new 
best-fit regions appear. The effect of channeling can easily be understood: 

• In CoGeNT, the channeled best-fit region has a DM mass reduced by a factor ^ 0.5 
(the square root of the quenching factor in eq. (3.5)). 

• In DAMA, scattering on lighter Na (on heavier Iodine) dominates if channeling is negli- 
gible (significant), and the best-fit region almost corresponds to the same DM mass, as 
shown in Fig. 13a. 

Adding channeling only in DAMA allows to get a best fit comparable to the un-channeled best 
fit. Fig. 13b shows one example of this possibility, assuming a 10% channeling in DAMA. 

Adding channeling also in CoGeNT moves its best fit to a smaller DM mass. Interest- 
ingly, with little channeling, no inelasticity and no isospin-violating interactions, we obtain an 
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channeling, 5 = -10.1 keV: 

CHpAjyA = 0.1 
CHcoGeNT = 

/„/,/;, = -1.48 
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----.. _?Na=a4 ; 




channeling 
CHqama = _ 
CHcoGeNT =0.1 

fpIL = 1 - 

Vfi = 220 km/s 
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?Na=0.4 : 
X^ = 105., 90.< 




Mdm in GeV 



Figure 13: Upper left: regions favored 
by DAMA (in magenta) and CoGeNT 
(in red) for (5 = and fp = fn shift 
when increasing the channeling fraction 
0% -^ 10% -^ 20%, as indicated by the 
arrows. Upper right; fit with inelastic 
dark matter and free fp/ fn os in Fig. 12, 
but adding a 10% channeling in DAMA. 
Bottom left; channeling in CoGeNT 
alone. As can be seen, 10% channeling 
allows for CoGeNT to evade all bounds 
while DAMA is still disfavored. Color cod- 
ing is described in Section ^.1. 



improved best fit (x^ ~ 100), corresponding to fitting CoGeNT compatibly with null experi- 
ments, but at the price of giving almost no signal in DAMA, which prefers a higher DM mass. 
We plot this possibility in Fig. 13c. 

4.10 Channeling plus inelasticity 

The incompatibility between the channeled DAMA and CoGeNT best fits in Fig. 13a or c, 
where we assumed fp/ fn = 1 and 6 = 0, can be eliminated varying these two parameters. Two 
kind of best fits appear: 

• The 'usual' best fit at higher Mdm ~ 6 GeV, with a slightly improved x^? as shown in 
Fig. 14a. 
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Figure 14: Fit with channeling in both DAMA and CoGeNT together with inelasticity. Left: 
global best fit at larger mass. Right; new best fit at smaller mass as discussed in Section 4-10. 
The dark magenta line is the bound from DAMA unmodulated signal. Bottom: DM predictions 
at the best fit point in the right panel. 

• A new kind of best fit, shown in Fig. 14b, that makes use of the channeled CoGeNT 
best fit at smaller Mdm ~ 3 GeV, and that employs inelasticity and channeling to also 
shift the DAMA best fit region to such a small mass. Indeed, a new local minimum of the 
DAMA-only fit appears at small mass for large enough channeling in DAMA. In Fig. 14b 
we assumed a 20% channeling fraction, that is enough to have both DAMA best fits with 
comparable x^- We find the exoDM (shown in the figure) to have a better global fit, 
however under the assumptions above, a lower mass region exists for endothermic DM 
too. 

The quality of the new best fit is now very significantly improved. It makes use of channeling 
only in DAMA and CoGeNT, of inelasticity and of isospin violation. It relies on channeled 
sodium events to constitute the bulk of the DAMA signal. 

Due to the small fractional modulation (the ratio between the modulated and unmodulated 
signal) characteristic of exothermic DM, non trivial bounds on this region come from demanding 
that the DM rate does not exceed the total rate measured by DAMA, in particular in the 1-2 



26 



keVee range. The dark magenta curve in Fig. 14b shows such bound. This constraint becomes 
less stringent for smaller values of \S\ (which worsen the fit) or under the assumption that the 
channeling fraction drops at energies below ~ 2 keVee. 

5 Conclusions 

We explored the prospects of dark matter for jointly explaining the modulations signals ob- 
served by DAMA and CoGeNT while complying with the constraints from other direct detec- 
tion experiments: XenonIO, XenonIOO, CDMS-Ge, CDMS-Si and SIMPLE. In the minimal 
scenario of spin-independent isospin-symmetric scattering there is a clear incompatibility, if we 
adhere to the data analyses performed by the experimental collaborations. Taking into account 
uncertainties on the astrophysical velocity distribution of DM and allowing for a higher sodium 
quenching factor in DAMA, gNa = 0.4 does not significantly improve the fits. Furthermore, 
DAMA and CoGeNT are mutually incompatible. 

Assuming different DM cross sections on proton and neutrons allows to significantly improve 
the situation, making DAMA and CoGeNT compatible and weakening one of the bounds 
from the null experiments. The best fit is for fp/fn ~ —1.5, which allows to weaken the 
dominant bounds from experiments with Xenon. However the overall global fit remains poor, 
see Fig. 8. Some additional improvement (although not very significant) is obtained by allowing 
for a DM/nucleon form factor that depends on the transferred momentum and/or on the 
relative velocity. A similar results holds for inelastic scattering of DM with a mass splitting, 
S = M^j^ — Mdm- Both up-scattering {6 > 0) and down-scattering {6 < 0) were considered. 
The best global fit to all relevant experiments is obtained as in Fig. 12 for \6\ ~ 10 keV. 

A modest channeling in DAMA allows a further slight improvement. Allowing also for a 
~ 10% channeling in CoGeNT results in a fit at smaller Mdm ~ 3GeV, that, together with a 
~ 20% channeling in DAMA and inelasticity, allows a good global fit (Fig. 14b). 

We stress the alternative simple possibility of mild channeling in CoGeNT alone: it allows 
for a good CoGeNT fit (without assuming a tuned inelasticity nor isospin-violating interac- 
tions) but excludes the DAMA signal. 

In conclusion, we find it hard to explain the full set of current direct detection experiments 
with spin-independent scattering of dark matter, and contemplate on the possibility that one 
or both of the positive signals does not arise from DM interactions. 

Note added: While this work was being finalized, related works appeared [62, 63]. 
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